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We present exact ground states in spin models with orbital degeneracy in one and higher 
dimensions. A method to obtain the exact ground states of the models when the Hamiltonians are 
composed of the products of two commutable operators is proposed. For the case of the spin-1/2 
model with two-fold degeneracy some exact ground states are given, such as the Valence-Bond 
(VB), the magnetically ordered, and the orbitally ordered states under particular parameter 
regimes. We also find the models with the higher spin and degeneracy which have the new types 
of VB ground states in the spin and the orbital sectors. 
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Quantum spin systems with orbital degeneracy 
have^a±lj-j|ctcji^mix_th£»retical and experimental intcr- 
estsMElM'ElBQl'BBy) These systems have rich phases 
such as the spin order, the orbital order, the spin liquid, 
and the orbital liquid. Examples of these systems are the 
Mott insulator of the transition metal (Mn, Cu, Cr, etc.) 
oxides,aB' the dynamical Jahn- Teller jjiolecular systems 
with a phonon coupling (TDAE-C6o)|3'tLand the heavy 
fermion systems in the insulator phaseJlJ' The strong 
interplay between the orbital and the magnetic ordering 
is interesting. Other interest is the existence of the spin 
and the orbital disorder states which are induced by the 
frustrations with the large quantum fluctuation arising 
from the spin-orbital interactions. There might already 
exist a physical realization of a auantum spin-orbital liq- 
uid in three dimension: LiNi02J3' The simplest model of 
the electron system with orbital degeneracy is the two- 
fold degenerate Hubbard Hamiltonian with Hund's rule 
coupling. In the strong coupling limit, the charge excita- 
tion has a large gap at 1/4- filling and this Hamiltonian 
is an effective model in the insulator phase such as the 
above spin systems. 

In the present letter, we consider the spin-orbital mod- 
els and give the exact ground states for the spin-1/2 and 
two-fold degenerate case. We extend the model to more 
than three-fold degenerate and half and higher spin cases 
which facilitate the higher dimensional VB states under 
certain parameter regimes. These VB states are the new 
types of the VB ground states and kinds of the Resonat- 
ing Valence Bond (RVB) states. 

First, we consider the Hamiltonian is defined by 



H = ^(W,,-A^')(V.,,-AS), 



(1) 



(^j) 



where the sum runs over all nearest-neighbor bonds 
(i,j). Hi J and Vij denote some Hermitian local op- 



by 1^) — |u) (g) |v), where |u(v)) denotes a state vector in 
u(v) subspace. The local Hamiltonian is positive semi- 
definite, and hence the total Hamiltonian is also positive 
semi-definite. The following two types of states are the 
possible ground states of the Hamiltonian. 

Type(i). — One possible ground state is written by 
\^g{A,B)) = |$"(A)) (g) |$^(B)). |$"(A)) and \^^{B)) 
are given by 






\^''{B))= n l*'^')o'(2) 

(»J)6-B 



where |«,j)" 



,(n),u(v). 



and 



/o - Y.nC'n'\<t>r)Z'{n = 1,...,7V) 
\va )ij is the lowest eigenstate of the local operator 
U{V) on sites i, j, whose degeneracy is N . A and B de- 
note a bond- covering (a configuration of the bonds) in 
u and V subspaces, respectively, which are constrained 
hy A\J B = { all bonds }. One site has z-^ and z^ 
bonds on the bond-covering A and S, respectively, obey- 
ing z < -^u + -Zv, where z denotes a lattice coordination 
number. Obviously, H|*g(^,-B)) = 0, thus |*g(^,-B)) 
is the ground state, because the eigenvalues of the Hamil- 
tonian are positive semi-definite and Y^q{A,B)) is the 
lowest eigenstate with an eigenvalue of zero. The ground- 
state degeneracy depends on the number of the bond- 
coverings. 

Type(ii). — If the state |$"(^') in one sector obeys 
H|$"(v)) ^ 0, where |$"(^)) = Hfe) I*, J>o^''\ the ground 
state of the global Hamiltonian is given by |^g) = 
|$u(v)^ g) |v(u)), where |v(u)) denotes an arbitrary state 
in v(u) space. Because Ti. is positive semi-definite and 
Tij^o) = 0, and hence |^g) is the lowest eigenstate 
with an eigenvalue of zero. In addition, the ground 
state of the Hamiltonian H ^ H + Q;7iv(u) (o^ ^ 0) 
is given by |*g) = 1$"'"^) «) |$^(")), where |$^(")) is 
the ground state of the 7iv(u) with the ground-state en- 

v(u) 



erators whose lowest eigenvalues are given by A^ and ergy ^^ '^^ since H - a^o is positive semi-definite and 



Q, respectively. 



If Uj i and V, , commute with each 



other: [Z^,;j, V^j] — 0, the eigenstate of H. is described 



""^j 



'«j 



{n ~ a^o'"^)|*G) = 0, and hence |*g) is the lowest 
eigenstate with the eigenvalue aE^ . 

Moreover, the ground states of the Hamiltonian: Ti — 
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E(,,,) an,ra{U^,, " KY^V.^j - AJ)" (a„,™ > 0, n, m > 1) 
are also the same as those of the two above cases. 

Using the above resuhs we discuss the ground states 
of the spin-orbital Hamiltonian: 

H = ^ JS^ ■ S, + Jrir^r^ + T^T^ + Arr^r/) 

+ J..(S, . S,)(rf r; + ryr] + ^t^t^), (3) 

where S^ denotes the spin operator and rf (a = x, y, z) 
denotes the a-component of the pseudo-spin operator. 
The pseudo-spin + and — (r/ — 1/2 and —1/2) cor- 
respond to the doubly degenerate orbital levels. This 
Hamiltonian is derived from the two-fold degenerate 
Hubbard model with Hund's rule coupling in the strong 
coupling limit at 1/4-filling by the second-order per- 
turbation.El' The Hamiltonian can be rescaled by \Jst\ 
and let us set | Js^-j = 1, lAi_j = Jsr^i ■ Sj, and Vi.j = 
TfTj -I- t^t}- + I^'t^Tj . The lowest eigenstate(s) of V^j is 
|singlet)[j — (|-|-—)ij- — | — -f)ij)/-\/2 with the eigenvalue 

— -^ — i for A' > —1, and are | + -f )ij and | )ij with 

the eigenvalue ^ for A' < — 1. For Jst — 1 the lowest 
eigenstate oiUi^j is singlet: |singlet)f .,- = (| ti)i,j - I it 
)ij)/-\/2, with the eigenvalue — |. We can find the fol- 
lowing exact ground states in one and higher dimensions. 
Dimerized ground states in one dimension. — For Jg = 
4^ -h 5 (A' = A > -1) and Jr = |, the Hamihonian 
(3) can be written as 7Y = H{i,j)^i,] + j){^i,j + ^ + 
2)+constant. The ground states are type (i) and given 

by 



where a ^ Js 



^' (A' < -1), /? = J^ - I > 0, and 



I*g) = n Isinglet)^_,+i ® |singlet)J_j-+i, 



(4) 



where i =even(odd) and j =odd(even). The ground 
states are the products of nearest-neighbor singlet dimers 
and two-fold degenerate. The translational symmetry is 
spontaneously broken although the interaction is only 
nearest-neighbor and no bond-alternation. The ground- 
state energy is Eq = — ^(A-(-2)L, where L denotes the 
number of sites. At A = 1, Kolezhuk and MikeskaO 
have already given the exact ground states in the lad- 
der model with a leg-leg biquadratic interaction by the 
matrix product ansatz. In the limit A — > cxd, the lowest 
states of the local Hamiltonian in the orbital part are 



two-fold degenerate (| 



/ij+i 



and I - +)j,j+i). The 



ground states of the global Hamiltonian in the orbital 
space are also two-fold degenerate and the orbital order 
(the Neel order in the pseudo-spin part) exists. At this 
point the system has a first-order phase transition. In 
other words, if an infinitesimal xy-element exists in the 
pseudo-spin space, then the ground states are the sin- 
glet, disordered states. The interaction of the transverse 
direction between spins and orbitals drastically enhances 
the quantum fluctuation. 

Ferromagnetic ground state in the pseudo-spin 
sector — When the Hamiltonian (3) can be expressed as 



n 



+ /5H^xz(7) 



3 A' 



aSi ■ Si 



7 < — 1 and Ti-xxzil) denotes the Hamiltonian of the 
XXZ model in the pseudo-spin space and is defined by 



Wxxz(7) - E(^^^/ 



rfrj+jr^r^), 



(5) 



1, and A 



^(A'-7)+7 



namely, for Jt > j, Jst - ^, "^^^ " - 37, 
(A' < —1), the ground state is type (ii) and is a fer- 
romagnetic state in the pseudo-spin space. In the spin 
space the ground state is equivalent to that of the anti- 
ferromagnetic (AF) Heisenberg model for Jg > — A'/4 
and of the ferromagnetic Heisenberg model for Jg < 
— A'/4. For the ferromagnetic spin states the ground- 
state energy is given hy Eq — —j [Jg + Jr7 -I- A' — I7] L 
in any dimension. In D = 1 the AF Heisenberg 
model is solved by the Bethe_Ansatz and has unique, 
massless singlet ground state.E^ The one-dimensional 
ground-state energy for Jg > — A'/4 is given by 
Eo = -i[^s(41n2-l) + J,7 + A'(ln2+i)-|7]L. 
For D > 2 the existence of an antiferromagnetic long 
range order(AFLRO) iru-the AF Heisenberg model has 
been rigorously proved,M but in Z? = 2 there is no 
proof of the existence of the AFLRO. Adding an in- 
teraction term S^,- -^ S^Sj {S > 0) to the Hamilto- 
nian at Jg = — J- and Jr — j, the ground states in 
the spin sector are equivalent to those of the AF Ising 
model. Thus, the system has the AFLRO on the bi- 
partite lattices in arbitrary dimensions with infinitesi- 
mal anisotropy of the z-direction. On the other hand, 
at Jg = — — and Jr = j the spin liquid states are 
induced by the infinitesimal perturbations such as an 
inhomogeneous interaction or a frustration reflecting a 
lattice topology. For example, let us consider the 1/5- 
dcplcted square lattice—which has the same topology 
as the CaV409 lattice,ll3 and consists of the four-site 
squares and the bridge bonds among the squares. If one 
adds an interaction term Ji I](ij)esquarcs ^« ' ^j ('^1 > 0) 
to the Hamiltonian, the exact ground state is a product 
of the four-spin singlet plaquette RVB (PRVB) states on 
the squares. On the other hand, adding an interaction 
term J2 1](ij)ebridgcs ^^ ' ^j ('^2 > 0) to the Hamilto- 
nian, the exact ground state is a product of the sin- 
glet dimer states on the bridges. We can exactly con- 
struct the PRVB and the dimer states without the next- 
nearest-neighbor interaction. An orbital long range or- 
der (ORLO) makes the systems have the spin liquid state 
under certain parameter regions. 

Ferromagnetic ground state in the spin sector — Next 
we set Jsr = — 1 and the lowest eigenstates of Uij are 
triplet whose eigenvalue is 
(3) can be expressed as 



J. When the Hamiltonian 






1, 



;m. 



A' 1 
-+2) + «S. 



+ PTi-xxzil) + constant, 
where a = Js + ^ + ^ < (A' > -l)and/3 = J^-i > 0, 



constant, 



2 — ^ \" — ^) """ A- "^r ^ 4 

namely, for Jg < -^ ~\, Jr> \, and A = ^^ -f 7 
(A' > —1), the ground state is type(ii) and a ferromag- 
netic state in the spin sector. In the orbital sector the 
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ground state is equivalent to that of the XXZ Hamiho- 
nian Ti-xxzil)- On a square lattice the AF XXZ model 
has the AFLRO for 7 > 1.66 that has been proved by 
Kubo and Kishipand-Nishiniori et.ai, using the infrared 
bounds method.ESElP This AFLRO corresponds to an 
OLRO in the orbital space. The state with the OLRO is 
what the electrons alternately occupy two orbital levels. 
Next, we consider the model with the spin-i and the 
pseudo-spin- 1 in one and higher dimensions described 

by 






+ /3{n-T,f + l]. (6) 



This model is regarded as the effective model with the 
spin-i for the 3-fold orbital degeneracy. For simplicity, 
we discuss only the isotropic case. 

In one dimension for a = 3/4, < /3 < 1/3, and 
7 = 2 — 4/3 the exact ground states are type(i) and 
the singlet dimer states in both the spin and the or- 
bital spaces, because one can write the Hamiltonian as 
^ = T,i^,J),n «»[S« • Sj + 3/4] [t, • Tj + 2]", where a„ > 
and n > 0. These states are two-fold degenerate and the 
translational symmetry is spontaneously broken. For the 
case of a = 3/4 and /3 = 1/3, the orbital ground state 
is equivalent to a Valence Bond Solid IVBS) state of 
the AfHeck-Kennedy-Lieb-Tasaki model.H In the spin 
sector the ground state is equivalent to that of the AF 
Heisenberg model for 7 > 2/3 and of the ferromagnetic 
Heisenberg model for 7 < 2/3 (type(ii)). At /3 = 1/3 
and 7 = 2/3 in the orbital sector the VBS state and the 
dimerized states are degenerate and in the spin sector 
all states are degenerate, and hence this point is multi- 
critical. 

At (a,/3, 7) — (|, |, |) the Hamiltonian is reduced to 

n = 2i:(,,) n^f 'Vir'^- Here T'ff^) and p(;=^) are 
the projection operators onto the subspace of the total 
spin S — 1 and the total pseudo-spin r = 2 on the two 
sites i,j, respectively. In arbitrary bipartite lattices with 
the coordination number three, e.g., a 2-leg ladder, a 
hexagonal lattice, a 1/5-depleted square lattice, we can 
find that the ground states are type(i) and given by 

|*g(A B)) = iDimer)" ® |VB)^, (7) 

|Dimer)^= J] i<^lb] ^ bla])\0), 

\YBy= n (44-44)10)' 

(^,f)eB 

where aj(6j) and cj(dj), denote the Schwinger boson's 
creation operators that create the S = 1/2 up(down) 
spin and the r = 1/2 -|-(— ) pseudo-spin at the site i, 
respectively. A and B denote a configuration of dimers 
(zs — 1) and a configuration of VB's (zr — 2) which has 
two bonds per site, respectively, constrained as A\J B = 
{all bonds} and A{^B — 0. The ground-state degeneracy 
is equal to the number of the dimer-coverings. On a lad- 
der and the higher dimensional lattices the number of the 
bond-coverings is infinity as the lattice size increasing to 
infinity. These ground states are kinds of the RVB states 
and the new types of the VB states induced by the in- 



terplay between the spin and orbital degrees of freedom. 
Taking some additional interaction terms to the Hamil- 
tonian, one can restrict the number of the degenerate 
states. The ground state of the Hamiltonian: H ^ H + 



^(ij)erungs y i,j 4 ^ ^W Z^(i j)elcgs y i,j 3^ 

x,y > is a product of all rung dimers in the spin sec- 
tor and two decoupled VBS states in the orbitaLsector. 
Therefore, the excitation energy has a finite gap.ES' An- 
other Hamiltonian: Ti. ^ H + xJ2m '^Moi^ ^ 0) where 
TC^Q denotes the Hamiltonian of_the Majumder- Ghosh 
model on the chain m{m = 1,2).ES' The ground states 
are four- fold degenerate and have two types. One is a 
checkerboard-type product of the singlet dimers along 
the legs in the spin sector and the VBS state connect- 
ing a chain in the orbital sector (fig. 1(a)). The other is 
a product of the two neighboring dimers along the legs 
in the spin sector and a product of the 4-site plaque- 
ttc VB states in the orbital sector (fig. 1(b)). These are 
the localized states as the dimer, the VB, and the finite 
VBS states. In the case of above two Hamiltonians, it is 
expected that each excitation energy has a finite gap. 

In general, the model can be easily extended to the 
higher spin and pseudo-spin cases. We can construct the 
models with the spin-5 and the pseudo-spin-T on the 
lattices whose coordination number is z, which have the 
following exact VB ground states: 

\^G{A,B)) = \<i>%A))<E>\^^{B)), (8) 

\^^{A))= n (44-4«])''^io)' 

ihj)eA 
\'^^iB)}= n i4d]-4clr^\Q). 

Here Ms and Mt obey Mg — 2S/zs — integer and M^ — 
2T/zr — integer, respectively, and z satisfies z < Zg + Zr- 
A{B) denotes a configuration of the bonds which has 
Zs{zr) bonds per site and Ayj B = { all bonds }. The 
Hamiltonian with these ground states is given by 

2S 2T 

^ = E E E kjM^-'^''^^^''' 

ii,j) J=Jo I=Io 



V^'.'"' (9) 



where Jq = 2S-Ms + l, h = 2T-i\/^ + l, and Kjj > 0. 
The simplest example of the one-dimensional models is 
given by H = E(,,,) [^^ ' ^3 + ^(^+1)] [n ■ Tj + T{T+1)] , 
whose ground states are the singlet dimer states with 
two-fold degeneracy in both the spin and the pseudo- 
spin spaces, where (z, Zs,Zt) = (2,1,1). Moreover, one 
can easily extend to the anisotropic cases. It is natural 
that the SU(2) symmetry is broken in the orbital sec- 
tor. The models have the same ground states of the 
isotropic cases as long as the lowest level of the two-site 
local Hamiltonian does not cross the higher levels. 

In summary, we have proposed a method to obtain the 
exact ground states of the models which are composed 
of the products of two commutable operators. Apply- 
ing this method to the quantum spin models with the 
orbital degeneracy as the Hamiltonian (3), we have ob- 
tained various ground states of such models and found 
that these models have the following interesting proper- 
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ties under particular parameter regimes. The frustration 
arising from the spin-orbital interaction of the transverse 
direction increases the quantum fluctuation and induces 
the singlet dimer ground states. The long range orders 
in one sector enhance various instabilities in the other 
sector in any dimension, e.g., the orbital orders make 
the systems have the spin liquid states and the AFLRO 
exist with infinitesimal anisotropy of the z-direction on 
the bipartite lattices. In addition, we have constructed 
models with arbitrary spin and orbital degeneracy which 
have the exact VB ground states whose degeneracy is in- 
finity in the thermodynamic limit on a ladder and higher 
dimensional lattices. These VB states include the corre- 
lations between the spin and orbital degrees of freedom 
and are the new types of the RVB states in D > 1. The 
exact ground states are obtained under restrict param- 
eter regions, but we expect that these states are adia- 
batically connected with the ground states in the large 
parameter space. 

The author would like to thank M. Takahashi, T. 
Kawarabayashi, M. Nakamura and N. Muramoto for 
their encouragement and helpful comments. 
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Fig. 1. Two types of the VB states on the ladder 
text). Pairs in the spin sector indicate dimer states. 
in the orbital sector denote (a) VBS, (b) finite VBS 



lattice (see 
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states. 



